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Abstract 

In this work a generalization of Onsager-Machlup's theory of time-dependent thermal fluctua- 
tions of equilibrium systems is proposed, to the case in which the system relaxes irreversibly along 
a non-equilibrium trajectory that can be approximated as a sequence of stationary states. This 
generalization is summarized by a canonical description of the dependence of the two-time corre- 
lation function C(i + T, t), and of the equal-time correlation function a(t) = C(t,t) (the covariance 
of the fluctuations), on the non-equilibrium relaxation time t 
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I. INTRODUCTION 



In this paper a generalization is proposed of some aspects of Onsager-Machlup's theory 
of equilibrium thermal fluctuations |2j, to conditions outside the so-called linear regime of 
irreversible thermodynamics 4], in which this theory is assumed to be universally valid 0, 
6fl . Although the original motivation of this work was to extend a theory of colloid dynamics 

on 

to describe the non-equilibrium slow dynamics in glass- forming colloidal systems [3, |8J, 
the resulting generalization actually constitutes a general canonical theory to describe the 
dependence of the two-time correlation function C(t + r, t) and of the equal-time correlation 
function o~(t) = C(t, t) (the covariance of the fluctuations) on the non-equilibrium relaxation 
time t. 

The dynamic properties of colloidal dispersions has been the subject of sustained interest 



for many years 
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These properties can be described in terms of the relaxation 



of the fluctuations 5n(r, t) of the local concentration n(r, t) of colloidal particles around its 
bulk equilibrium value n. The average decay of 5n(r,t) is described by the two-time cor- 
relation function F t (k,r) = (5n(k,t + t)5ti(— k, t)) of the Fourier transform <5n(k, t) of the 
fluctuations 5n(r,t), whose equal-time limit is S t (k) = F t (k,r = 0) = (5n(k, t)Sn(— k, t)). 
One may refer to the time r as the correlation time, and to the time t as the relaxation 
time. If some external (or internal) constraints that kept a system at a certain macroscopic 
state are broken at time t = 0, the system relaxes spontaneously to its new thermodynamic 
equilibrium state, and one may also refer to t as the waiting or ageing time @, S| • If the 
system has fully relaxed to, and/or remains in, a thermodynamic equilibrium state, these 
properties no longer depend on t, i.e., F t (k,r) = F(k,r) and S t (k) = S(k). The equilibrium 
stationary correlation function F(k, r) is referred to as the intermediate scattering function, 
and its initial value S(k) as the static structure factor. These properties can be measured 
by a variety of experimental techniques, including (static and/or dynamic) light scattering 

3, B,Q- 

The static structure factor S(k), being a thermodynamic property is amenable to theo- 
retical calculation using statistical thermodynamic methods [12J. The fundamental under- 
standing of F(k,r), on the other hand, requires the development of theoretical methods to 
describe the diffusive relaxation of the local concentration fluctuations, and a number of such 
approaches have been proposed for their theoretical calculation 
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One of the earliest was developed by Hess and Klein 
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mode-coupling self-consistent theory of molecular liquids 



13 1, who translated to colloids the 
18, Q. More recently, Nagele and 



collaborators elaborated further this mode-coupling theory of colloid dynamics 
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An independent alternative theory of the dynamic properties of colloidal dispersions has 
been developed within the last decade, and is referred to as the self-consistent generalized 
Langevin equation (SCGLE) theory of colloid dynamics Q, H Q, S, Q- This theory 
has been recently applied to the description of dynamic arrest phenomena in several spe- 
cific colloidal systems that include mono-disperse suspensions with hard-sphere interactions, 
moderately soft-sphere and electrostatic repulsions, short-ranged attractive interactions, and 
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35], the 



model mixtures of neutral and charged particles 

Until recently, and in spite of the long tradition in the study of glasses 
only well-established and successful theoretical framework leading to quantitative predictions 
of the glass transition was the conventional mode coupling theory (MCT) of the ideal glass 



transition 
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35 



36]. Many of the predictions of this theory had been systematically 



confirmed by their detai 



systems [37 
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ed com par ison with experimental measurements in model colloidal 
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46| . In this context, we can mention that the 



more recently-deve 
scenarios as MCT 



oped SCGLE theory of dynamic arrest leads to similar dynamic arrest 
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26j for several specific (mostly mono-disperse) systems, although for 



colloidal mixtures important differences may appear in some circumstances, as reported in 
Refs. 



30 



31 



32|. 



An important common feature of both theories in their current status is that they are able 
to predict the regions of the control parameter space in which the system is expected to be a 
glass, i.e., it predicts what we refer to as the "dynamic arrest phase diagram" of the system 
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32j . Since the predicted state only depends on the end value of the control parameters, 



these predictions only apply to "reversible" or "equilibrated" glasses. Thus, while it is im- 
portant to pursue the application of these two theories to specific idealized or experimental 
model systems and to compare their predictions, it is also important to attempt their ex- 
tension to the description of "non-equilibrated" glasses, for which no "dynamic arrest phase 
diagram" will make sense without the specification of the detailed non-equilibrium process 
leading to the apparent end state. Ageing effects, for example, should be a fundamental 
aspect of the experimental and theoretical characterization of these non-equilibrium states. 
These preoccupations have been addressed in the field of spin glasses, where a mean-field 
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theory has been developed within the last two decades 47J. The models involved, however, 



lack a geometric structure and hence cannot describe the spatial evolution of realglass for- 
mers. Thus, although computer simulations 48 



491 ] and experimental studies 
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have 



provided important information about general properties of ageing, no quantitative theory 
is available so far to describe the irreversible formation of structural glasses. 

Almost a decade ago, however, an attempt was made by Latz [50] to extend MCT to 
describe the irreversible relaxation, including ageing processes, of a glass forming system af- 
ter suddenly driving it into the glassy region of its dynamic arrest phase diagram. A major 
aspect of this work involved the generalization to non-equilibrium conditions of the conven- 
tional equilibrium projection operator approach [5JJ to derive the corresponding memory 
function equations in which the mode coupling approximations could be introduced. 

In very recent work [52] a similar extension of the SCGLE theory of dynamic arrest has 
been proposed by the present author. In contrast with MCT, the fundamental basis of the 
SCGLE theory does not involve the use of projection operator techniques. Instead, it is 
based on Onsager-Machlup's theory of time-dependent thermal fluctuations of equilibrium 
systems. Thus, its extension to non-equilibrium conditions calls for a generalization to the 
nonlinear regime of the general and fundamental laws of linear irreversible thermodynamics 
and the corresponding stochastic theory of fluctuations, as stated by Onsager {3], 4] and by 
Onsager and Machlup [1, 21, respectively, with an adequate extension 53, 54 1 to allows the 



description of relaxation phenomena involving memory effects. In Ref. 52j this general- 
ization of Onsager's theory was only outlined. The main purpose of the present paper is 
then to explain in more detail the main features and the underlying assumptions of such 
generalized canonical theory of time-dependent fluctuations around the irreversibly relaxing 
state of non-equilibrium systems. 

For this, and since the literature on non-equilibrium extensions of equilibrium theories 



is quite diverse and extensive , [55|, and in order to normalize the basic concepts and 
scientific context, we emphasize that we follow to a large extent the philosophical approach of 
Joel Keizer's statistical thermodynamic theory of non-equilibrium processes [5], particularly 
his account of the Onsager picture (chapters 1 and 2 of Ref. 5j) and some aspects of 
Keizer's extension of this picture to non-equilibrium (chapters 3 and 4 of the same reference). 
However, for our present purpose we do not find necessary to adhere to Keizer's detailed 
mechanistic statistical description in terms of elementary processes. 
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Thus, taking some elements from Keizer's theory and other elements not considered by 
him (such as the non-Markovian extension of Ornstein-Uhlenbeck processes 53|), the present 
paper proposes a generalization to nonlinear conditions of some elements of Onsager's canon- 
ical theory. We refer essentially to the time-evolution equations of the first two moments of 
the conditional probability distribution (i.e, the mean value and the covariance of the state 
variables) and of the two-time correlation function. These time-evolution equations will be 
derived directly from the postulated linear stochastic equation with additive noise, rather 
than from the time-evolution equation of the conditional probability distribution, which is 
another possible route [561 ]. 

Onsager-Machlup's theory contains elements of purely mathematical nature that must be 
clearly distinguished from the purely physical assumptions about the behavior of physical 
systems. Among the former, the most relevant is the concept of Ornstein-Uhlenbeck pro- 
cess, i.e., a stationary Gaussian Markov stochastic process 



57 



. This is the mathematical 



model employed to state the main assumption of physical nature in Onsager-Machlup's the- 
ory, which can then be stated by saying that the spontaneous fluctuations of an equilibrium 
system can be described as an Ornstein-Uhlenbeck process. Thus, after a brief qualitative 
discussion in section |TT] on the assumed general features of the thermodynamic conceptual 
framework that we have in mind to describe the evolution of the system, section II III summa- 
rizes the definition and the main properties of an Ornstein-Uhlenbeck process as a mathe- 
matical model. In Section HVl we employ these concepts to summarize Onsager's theory (with 
this term we refer to Onsager-Machlup's theory plus Onsager's own work on the properties 
of the linear relaxation equations of irreversible thermodynamics (si, 4\). The generality of 
the underlying mathematical model makes Onsager's theory applicable to conditions not 
considered in the original work. These extensions serve as the basis for its generalization to 
the non-linear regime, which we present in Section |VJ To simplify the reading of this paper, 
in Section I VI I we summarize the main concepts and results that constitute this generalized 
canonical theory and in Section IVHI we provide an illustration of its concrete use by re- 
viewing its application to colloid dynamics. The significance of the results of this particular 
application of the generalized canonical theory is also discussed in the final section IVIIII 
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II. NON-EQUILIBRIUM RELAXATION AND ENTROPY LANDSCAPE 



Let us consider a system whose macroscopic state is described by a set of C extensive 
variables Oj(t), i = 1, 2, C, which we group as the components of a C-component (column) 
vector a(i). The fundamental postulate of the present statistical thermodynamic theory 
of non-equilibrium processes is that the dynamics of the state vector a(£) constitutes a 
multivariate stochastic process, described by a deterministic equation for its mean value a(t) 
and by a linear stochastic equation with additive noise for the fluctuations 5&(t) = a(t) — 
a(£). The basic assumption is that the mean value a(£) coincides with the macroscopically 
measured value, and that its time evolution is described by an equation of the general form 

^j® = *[!(*)], (2-1) 

where the (generally non-linear and spatially non-local) functional dependence of the C- 
component vector 1Z [a(t)] on a(i) includes both, dissipative and mechanical (i.e., conserva- 
tive) terms j^. 

Under some conditions, the relaxation equation above may admit time-independent, or 
stationary, solutions denoted by & ss . This means that a ss solves the equation 

da ss 

— = ^[a* s ]=0. (2.2) 

Stationary states may result from the non-linear mathematical structure of the relaxation 
equation, Eq. (12. ip . These are "true" non-equilibrium stationary states, in the sense that 
their stability requires the continuous input and output of energy, matter, etc., thus involving 
intrinsically open systems. In our present discussion we shall not refer to this kind of 
stationary states. Instead, we refer to stationary states that correspond to absolute or local 
maxima of the total entropy, implying conditions for the thermodynamic equilibrium of the 
system within the constraints imposed by its interactions with external fields and reservoirs 
and within the internal constraints resulting from, and sustained by, strong intermolecular 
interactions, such as in meta-stable or in arrested states. 

For a closed system with fixed external constraints (i.e., with given time- independent 
external fields, including those of the isolating and confining walls), there is one state that 
corresponds to the absolute maximum of the function entropy S = S[a\ within the set T 
of macroscopic states a consistent with those isolation conditions. According to the second 




FIG. 1: Schematic illustration of the "entropy" landscape of an isolated system whose fundamental 
thermodynamic relation S = S[a] involves a single macroscopic variable a. The isolation condition 
constrains the possible states a to a subset T of the entire state space of the system; in the figure, 
this isolation domain T, is represented by the interval a m i n < a < a max . The absolute entropy 
maximum is located at a = a eq . Fig. (a) illustrates the extreme case in which this global entropy 
maximum is the only maximum, whereas Fig. (b) illustrates the case in which there are many 
additional local entropy maxima, one of them indicated by its location a = a ss 

law of thermodynamics, such state, denoted by & eq , is the thermodynamically most stable 
equilibrium state. Fig. 1(a) illustrates schematically the simplest scenario in terms of an 
"entropy landscape" that exhibits a single and "unquestionable" entropy maximum. This 
global maximum a 69 must then be the only attractive stationary state of the non-linear 
relaxation equation ( 12.11) . whose basin of attraction must then coincide with the whole 
subspace T of states that are consistent with the isolation condition. This means that if 
the system is initially prepared in any state a e T, its mean value a°(t) will evolve in time 
according to Eq. (12. ip relaxing eventually, but surely, to the equilibrium state a eq . 

Severe internal constraints may, however, give rise to very many additional equilibrium 
states, denoted generically as a ss , that correspond to local entropy maxima, as illustrated 
schematically in Fig. 1(b). In contrast with the previous scenario, in this case, depending 
on the location of the initial state a e T, the system will not necessarily evolve towards the 
absolute thermodynamic equilibrium state & eq , since there are now many competing local 
entropy maxima which also act locally (i.e., within its own basin of attraction) as attractive 
stationary states representing meta-stable thermodynamic equilibrium states. In fact, each 
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of these meta-stable states, such as that represented by a ss in Fig. 1(b), will actually be 
as stable as the "most stable" equilibrium state a eq , as long as no external (or internal!) 
perturbations drive the system, even momentarily, out of the basin of attraction T(a ss ) of 
this particular stationary state. We assume that the basins of attraction of each pair of 
different stationary states in T are disjoint, i.e., that T(a^ s ) f]T(al s ) = for any pair a\ s 
and a| s , and that the union of the basins of attraction of all these stationary states, including 
a. eq , equals the state space T of the isolated system, [J aSS T(a ss ) = T. 

In the previous discussion we have assumed that the entropy landscape in Fig. 1(b) 
remains constant in time so that if the system was prepared in the basin of attraction 
of a given stationary state, it will be trapped in this basin, eventually relaxing to the 
corresponding stationary state sl ss , in which it will remain indefinitely afterwards. This 
means that all the solutions of the relaxation equation in Eq. (14.61) consist of trajectories 
confined to the basin of attraction of one specific stationary state. In reality, it is impossible 
that such situation can be sustained indefinitely, since there are many possible mechanisms 
that allow the trajectory a°(£) to explore states outside the basin of attraction in which the 
system was initially prepared. If such excursion occurs to the basin of a stationary state 
with a still larger local entropy maximum, the second law of thermodynamics dictates that 
the system will proceed to the stationary state of the new basin of attraction. If this process 
occurs repeatedly, the system will eventually reach the most stable equilibrium state a eq , 
thus restoring to some extent the irreversible process described in the simpler conditions 
illustrated in Fig. 1(a), but with a probably much slower dynamics because the transitions 
from basin to basin may involve activation barriers that must await for the occurrence of 
adequate spontaneous fluctuations. In this process, the lifetime of these instantly stable 
stationary states may exceed the experimental timescales, thus giving rise to meta-stability 
conditions, and/or to conditions in which the system may appear to be dynamically arrested. 

A fundamental assumption of the present theory will be that the actual relaxation of a 
system under the more complex conditions illustrated in Fig. 1(b) can be approximated 
by a sequence a°(£ a ) = a s J (a = 0, 1, 2,...) of locally stationary states a^ generated by 
the recurrence relation a* s +1 = a^ + 1Z [aj s ] (£ Q+ i — £ a ), starting with a^ = a . Although 
these locally stationary states a ss are not stable equilibrium states, we will assume that 
their properties may be described by quite similar fundamental principles as the most stable 
equilibrium state a eq . Clearly, the conceptual framework and the approximations employed 
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in the theoretical modeling of the dynamics of a°(t) must reflect the fundamental difference 
between the two qualitatively different scenarios illustrated in Fig. 1. Onsager's linear 
irreversible thermodynamic theory of fluctuations has in mind the simpler scenario 

illustrated by the entropy landscape of Fig. 1(a). We are interested, however, in generalizing 
Onsager's theory to the conditions represented schematically by the more complex scenario in 
Fig. 1(b). The fundamental postulate of Onsager's theory is that the macroscopic dynamics 
of the system is not described by a deterministic equation for the state vector a(t). Instead, 
it is postulated that the macroscopic state of the system is described by a statistical physical 
ensemble whose mathematical representation involves the assumption that a(t) constitutes 
a multivariate stochastic process, and more specifically, a stochastic process referred to as an 
Ornstein-Uhlenbeck process. Thus, before explaining the main physical assumptions made 
in Onsager's theory, let us review the purely mathematical framework of stochastic processes, 
in which such physical assumptions can be stated most economically. 

III. ORNSTEIN-UHLENBECK PROCESS, A REFERENCE MATHEMATICAL 
MODEL 



15 
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57 



^his^ section we summarize the main concepts that define an Ornstein-Uhlenbeck process 
. As a mathematical object, a stochastic process a(t) is defined in terms of the 



joint probability density W m (a.i, ti, a 2 , t 2 ; a m , t m ) for the state vector a(t) to have a value 
in the interval a^ < a(tj) < a^ + dan for i = 1, 2, m. We say that this stochastic process is 
fully determined if we know the probability densities for all possible positive integer values 
of m and all possible sets of times (ti, t2, t m ). If the stochastic process is Markovian, 
however, a great simplification occurs, since in this case all these probability densities can 
be written in terms of only W 2 (ax, t\] a 2 , t 2 ). This probability density can be written as 
iy 2 (ai, t\\ a 2 , t 2 ) = PFi(ai, ti)P 2 (ai, ti | a 2 ,t 2 ) where P 2 (ai,ti | a 2 ,t 2 ) is the conditional 
probability density that a(t 2 ) has a value in the interval a 2 < a(t 2 ) < a 2 + c?a 2 provided that 
for sure a(ti) = a x . 

A stochastic process is said to be stationary if all its probability densities are time- 
translational invariant, i.e., if for all real values of s we have that W m (ai + s,ti;a 2 ,t 2 + 
s; a m , t m + s) = W m (sn, £i; a 2 , t 2 ; a m , t m ). Thus, if in addition to being Markovian the 
stochastic process is also stationary, then Wi(a.x,ti) = W(a.i) and P 2 (a 1 ,t 1 | a 2 ,t 2 ) = P(ai | 
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a2,^2 ~ ti). Assuming that lim 4 ^ 00 P(ai | a2,t) = W(&2), then we have that a stationary 
Markov process is determined solely by the conditional probability density P(ai | a2, t^ — 

The knowledge of this probability density is equivalent to the knowledge of all its mo- 
ments. A final great simplification occurs when we assume that the stationary Markov 
process is, additionally, Gaussian, i.e., such that all the moments of P(a | a, t) can be 
written in terms of only its two lowest-order moments as 

P(a° | a,t) = [(27r) c det ( r(t)]e-[( a - s0 W) toCT " 1 W o ( a - s0 ( i »]/ 2 , (3.1) 

where the dagger means transpose, the circle "o" means matrix product, and where the 
conditional mean value a°(£) and the covariance matrix a(t) of the fluctuations 5&(t) = 
a(i) — a°(£) are defined, respectively, as 

a°(t) = J aP(a | a,t)da (3.2) 

and 

a(t) = 6a{t)6rf{t) = J - a°(t))(a - a°(t)) t P(a | a, t)da. (3.3) 

We notice that in the long-time limit, P(ao | a, t) attains its stationary value VT(a) = 
[(27r) c deta ss ]e-[( a - s3S ) toff3 ^ lo ( a - sss )]/ 2 , with a ss = a°(t -> oo) and a ss = a(t -> oo). 

A stochastic process that is stationary, Gaussian, and Markov, like the one just discussed, 



has a proper name: it is referred to as an Ornstein-Uhlenbeck process |57l . |58| . From a 
mathematical point of view this is a most attractive and useful mathematical model, since 
its full determination is reduced to the determination of two objects, namely, the mean 
value a°(t) and the covariance cr(t), whose determination then becomes a central issue. This 
issue, however, is partially settled when one demonstrates that each of these two properties 
satisfies a deterministic relaxation equation whose structure is also determined by the model 
itself, as we now discuss. 



A. Relaxation equations for a°(i) and <j(t) 



Besides the definition given above, the Ornstein-Uhlenbeck 

n 

other alternative but mathematically equivalent definitions [5 



stochastic process can be given 
The most relevant for our 
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present purpose defines the stochastic process a(t) as the stationary ensemble constituted by 
all the possible "realizations" of the stochastic process, which correspond to all the possible 
solutions of a linear stochastic differential equation with additive noise for the fluctuations 
5a(t) = a(t) — ai ss . Such stochastic equation has the following general structure 

^ = W o*a(f)+f(f), (3.4) 

with TC being aCxC relaxation matrix and with the C-component vector f(i) being a "white 
noise", i.e., a stationary and Gaussian stochastic process which is, however, not Markovian, 
but ^-correlated, i.e., such that < f(t)f(t') >= 72<5(t — t') (with "<•••>" indicating the 
average over all the possible realizations of the noise f(£)). In addition, it is assumed that 
< f(£) >= and < f(t)Saf(0) >= 0. The stochastic process Sa(t) thus defined can be shown 
to be Gaussian, Markov, and stationary, i.e., to be an Ornstein-Uhlenbeck process. 

The main reason one might prefer to define a stochastic process in terms of a linear 
stochastic equation with additive noise, is that one can derive the deterministic time evolu- 
tion equations for a°(£) and a{t) much more directly from such stochastic equation than from 
the definition of these properties in Eqs. (13.21) and (13. 3p . which require the previous deriva- 
tion of the time-evolution equation of the full probability distribution function P(ao | a, t). 
To illustrate this, let us notice that the solution of Eq. (13 .4p can be written as 

5a(t) = e nt o 5a(0) + / e n ^ o f(r)dr (3.5) 



with <5a(0) = a — a ss . Notice also that the mean value 5&(t) = a (t) — a ss is the average 



of 5&(t) over the realizations f(t) but not over the initial conditions <5a(0), i.e., Sa(t) =< 



5&(t) > with the initial condition 5a(0) = a — a ss . Thus, taking the average of Eq. ( 13. 51) 
over the realizations f(i), it follows that the a°(t) is given by 

a°(t) =a ss + e M o [a - a ss ]. (3.6) 

As long as the real part of all the eigenvalues of the relaxation matrix 7i remain negative, 
this expression interpolates a°(t) between its arbitrary initial condition a and its asymp- 
totic stationary value a ss . Clearly, this expression for a°(t) is the solution of the ordinary 
differential equation 
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dAa.(t) . . . . , 

-^p^oAafi), (3.7) 

where Aa(t) = a°(i) — a ss , with initial condition Aa(0) = a — a ss . Thus, we conclude 
that in the Ornstein-Uhlenbeck process the stochastic equation for the fluctuations 5a.(t) 
coincides with the equation that governs the decay of the deviation Aa(t) of the conditional 
mean value a°(t) except for the additive white noise term f(t) in Eq. (13.41) . We also notice 
that the previous relaxation equation for a°(t), as well as its explicit solution in Eq. (13.61) . 
only depend on the relaxation matrix 7i for given initial condition a . 

The relaxation equation for the covariance a(t) can be derived in a similar manner. Thus, 
by multiplying Eq. (13. 5p on the right by its transpose, and averaging the resulting expression 
both, over the realizations of i(t) and over the initial values £a(0), it is not difficult to derive 
the following expression for cr(t), 

a(t) = e m o[a° + 2 f e~ UT o 7 o e'^dr} o e nh , (3.8) 
Jo 

where cr° is the given initial value of cr{t). By differentiating this expression with respect to 
time, er(t) is seen to satisfy the differential equation 

^ = ^0^ + ^)0^ + 27 (3.9) 
with initial condition <j(0) = a . 



B. The fluctuation-dissipation theorem 

The formal solution for a(t) in Eq. (13. 8p can be written as a(t) = e m oa°oe ntt + 2 J Q * e Hs o 
7 o e H]s ds which, in the limit t — > 00, reads a ss = 2 J °° e Hs 070 e nts ds. Subtracting this 
expression from the previous one, and using the change of variable r = s — t, Eq. ( 13.81) may 
be given the following alternative expression 

a(t) = a ss + e m o [a - a ss ] o e nh , (3.10) 

which interpolates the covariance a(t) between its given initial value a and its asymptotic 
stationary value a ss . Taking the time derivative of this expression and comparing it with 
Eq. ( I3.9p . finally leads to the so-called fluctuation-dissipation theorem, 
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Hoff ss + a ss oH t + 27 = 0. (3.11) 

This relationship between the matrices 7i, a ss , and 7 constitutes a necessary and sufficient 
condition for the process to be stationary, and also follows directly from Eq. (13.91) provided 
that the long-time stationary solution a ss exists. Besides guaranteeing stationarity, this 
relationship also sets stringent conditions, of purely mathematical nature, on the structure 
of the relaxation matrix TC. To see this let us define the matrix £ as the product 

C = -Hoa ss . (3.12) 

Then, the fluctuation-dissipation relationship in Eq. (13. lip can be written as 

j = ^Y^ = C t . (3.13) 

This equation can be viewed as a condition that determines the matrix 7 (a measure of the 
time-correlation function of the additive white noise f(i)) in terms of the symmetric part £ s 
of the matrix £ defined, in its turn, in Eq. (13.121) in terms of the relaxation matrix H, and 
the stationary covariance a ss . 

The same equation can also be viewed as a condition of the general structure of the 
relaxation matrix Ti, which then must be such that it can be factorized as 



H = -£oa ss -\ (3.14) 

with the matrix £ being identical to the symmetric matrix 7 plus, at most, an antisymmetric 
matrix £ a = {£ — C')/2. Thus, we conclude that the fluctuation-dissipation relationship 
dictates that the relaxation matrix 7i is not completely arbitrary, but must comply with 
the rather rigid format of Eq. (13.141) . For example, the most general structure of the linear 
relaxation equation for the mean value Aa(t) in Eq. ( 13. 71) must be of the form 

w = -rotr^oAaft), (3.15) 

a format which, in the case that the stationary state sl ss is the thermodynamic equilibrium 
state a. eq , will soon be identified with the canonical form of the relaxation equations of linear 
irreversible thermodynamics. 
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C. Time- dependent correlation function C{t',t) 



Another important property of the stochastic process a(t) is its two-time correlation 



matrix C{t' } t) = Sa(t')6a^(t) where the over-line refers to an average over all the possible 
realizations of the stationary noise f(t) plus the average over the distribution of the initial 
values <5a(0). We may evaluate this property by using again Eq. (13.51) evaluated at t' and 
at t to form the product 5a.(t')5a) (t) . Taking the average of this product, we may use the 



fact that f(t')tf(t) = 72<J(t' - t) and f(t)5at(0) = to show that C{t',t) may be written, 
assuming that t' > t, as 

C(t',t) =e Hl - t '- t) oa{t), (3.16) 
where the expression for a(t) in Eq. (13.81) was employed. 



Defining t' = t + r, one can also write C(t',t) as C(t +r,t) = 8a.(t + r)5at(i) = Ct(r) 
to describe the decay of the correlations with the correlation time r when the system has 
evolved during a relaxation time t after it was prepared in an initial state represented by 



the distribution of initial values 5a(0) with mean value 5a(0) = a — a. ss and covariance <t°. 
One can also derive the time evolution of C t (r) by writing t + r as the argument of Sa(t) in 
Eq. ( 13. 51) . whose derivative with respect to r for fixed t can be shown to be given by 

d6a ( t + T ^ = n o 5a(t + r) + f(t + r). (3.17) 
dr 

Multiplying this result on the right by the transpose of Sa.(t) and averaging over the real- 
izations of f(t + r), one is led to the following relaxation equation for C(t + r,t), 

dC(t + r,t) 
dr 

with initial condition C(t + 0, t) = o~(t). The solution of this equation is again given by Eq. 
( 13.161) . now written as 

C(t + r,t) = e Hr oa(t) (forr>0). (3.19) 

Let us notice that for sufficiently long relaxation times t after the initial condition set 
at t — 0, the process recovers its stationary condition, in which a(t — > oo) = a ss and 
C(t + T,t) = C(t, 0) = C ss {t), with 



C ss (T)=\imC(t + r,t) = e nT oa ss (for r > 0). (3.20) 

t— *oo 
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In summary, the Ornstein-Uhlenbeck stochastic process is fully determined by the time 
evolution of the mean value a°(t) and the covariance cr(t), governed by the deterministic 
equations (13. TR and ( 13.91) . The solution of these equations, given in Eqs. (13.61) and (13.101) . 
is a simple exponential interpolation between the given initial state (a , a ) and the final 
stationary state (a ss , a ss ) with relaxation times given by the eigenvalues of the relaxation 
matrix TC. This matrix, however, is not completely arbitrary, but must conform to the 
format in Eq. (13. 14[) . Nevertheless, this matrix becomes the kernel of the model, since it 
is the main ingredient in terms of which all the properties of the process are determined. 
This includes, besides a°(t) and u(t), the time- dependent correlation function C(t',t) just 
discussed, and written explicitly in terms of H in Eq. (13. or (13.191) . 



IV. ONSAGER'S THEORY OF EQUILIBRIUM THERMAL FLUCTUATIONS 

The mathematical infrastructure just reviewed can now be employed to described specific 
physical phenomena in a very efficient manner. This approach was inaugurated by Langevin 



59] with his celebrated equation for the thermal fluctuations of the velocity V(t) of a 



Brownian particle. In reality his work triggered the development of the mathematical field 



of stochastic processes [60| , from which we borrowed the concepts summarized in the previous 
section. Today, however, we may state his theory most efficiently as the assumption that the 
cartesian vector V(t) constitutes an Ornstein-Uhlenbeck process specified by the diagonal 
relaxation matrix TCij = —((/M)5ij (with i, j = 1, 2, 3, and with ( and M being the friction 
coefficient and the mass of the Brownian particle). From this assumption, all the results of 
the previous section apply with a(t) replaced by V(t), and with this specific value of 7i. 

With this identification of the relaxation matrix 7i, the textbook presentations of 
Langevin's theory of Brownian motion [12j then become just a review of the physical mean- 
ing, worded in the specific context of Brownian motion, of the mathematical properties of 
the underlying Ornstein-Uhlenbeck mathematical model summarized in the previous sec- 



tion. For example, according to Eq. (13.41) . and as Langevin puts it in his original work [591 ]. 
the thermal fluctuations V(t) are then governed by the stochastic version of this linearized 
Newton's equation, namely, M(d~V(i) / di) = — £V(t) + f(t). The random force f(t) is then 
assumed to be "indifferently positive and negative and with a magnitude such that it main- 
tains the agitation of the particle, which the viscous resistance would stop without it" , a 
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physically intuitive manner employed by Langevin to express his assumption that f(t) is a 
stationary delta correlated noise. 

It is then instructive to recall Langevin's arguments leading to the identification of the 
relaxation matrix TC above. These are extremely simple: the conditional mean velocity V(t) 
of a colloidal particle obeys the same equation that governs the motion of a macroscopic 
particle under the influence of the frictional resistance R[V(t)] of the solvent, i.e., Newton's 
equation, M(dV(t)/dt) = R[V(i)]. The functional dependence of the resistance R[V(i)] 
on the velocity V(i), referred to as the "constitutive" relation, is in general non-linear and 
poorly understood, but its linear approximation R[V(t)] ~ —(V(t) leads to the simplest 
form of the relaxation equation for V(t), namely, M(dV(t)/dt) = — £V(i). Comparing with 
Eq. (13. 7p . one immediately identifies the relaxation matrix to be Tiij = —((/M)Sij. 



A. Linear laws and linear regime of irreversible thermodynamics 

Onsager's theory of equilibrium thermal fluctuations is, to a large extent, an extension 
of Langevin's theory, in which the vector V(£) is extended to be a generic state vector 
a(t). Just like in Langevin's theory, the first main postulate of Onsager's theory is that the 
macroscopically observable state variables correspond to the mean value a(£) of a stochastic 
process whose time evolution is governed by a nonlinear relaxation equation that we write 
in the general form of Eq. (12.11) , namely, 

^&=n\m- (4-i) 

This general non-linear equation is assumed to represent, for example, the hydrodynamic 
transport equations or Fick's diffusion equation, as particular cases 0, |(5lj]. The generally 
nonlinear dependence of the "flux" IZi [a(i)] (the time rate of change of the extensive variable 
~a~i(t)) on the state vector a(i) is referred to as the constitutive relation, a relation that turns 
Eq. (14. ip into a closed equation for a(t). 

Although the explicit form of this relation is not known in general, in the so-called lin- 
ear regime, discussed below, the resulting relaxation equation must conform exactly to a 
rather strict format. According to the point of view adopted in the conventional account of 
Onsager's theory , |6l| , the format of the relaxation equations in the linear regime is 



a direct consequence of one of the most fundamental physical principles of non-equilibrium 
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irreversible thermodynamics, namely, the phenomenological so-called linear laws |6l|]. This 
principle states that the "flux" IZi [a(t)] is proportional to the C-component vector of "ther- 
modynamic forces", AF(t) = F[a(t)] — F eq , whose components describe the instantaneous 
imbalance of the conjugate intensive variables Fj[a(t)] = k^ 1 (dS[a]/daj) a= -^ with respect 
to their equilibrium value F eq . Thus, the linear laws state that 

K [a(t)] = -C* [a(t)] o AF(t) (4.2) 

with the proportionality matrix C* [a(i)] being referred to as the matrix of Onsager's "kinetic 
coefficients". We assume that this matrix is known in principle from phenomenological 
considerations, as in Fourier's heat conduction or in Fick's diffusion laws. The star on 
£* [a] is as reminder of the phenomenological nature of this property and of its dependence 
on the state vector a. Clearly, this principle requires the existence of the state function 
entropy S = S[a] which, for given conditions of isolation and given fixed external fields on 
the system, has its maximum at a particular state a eq , as prescribed by the second law of 
thermodynamics. Let us stress that the entropy S refers to the entropy of the closed system, 
in which case F[ q = 0; if the system is in contact with thermal and/or chemical reservoirs, 
the corresponding F[ q may, however, differ from zero. 

The linear laws described by Eq. (14.21) do not necessarily imply that the resulting re- 
laxation equation for the macroscopic variables a°(£) will in general be linear. Onsager's 
theory, however, is restricted to the so-called linear regime around the thermodynamic equi- 
librium state a. eq . Thus, imagine that the system is prepared in an initial state a that lies 
in the basin of attraction T(a eq ) of the equilibrium state. Then, after some time the tra- 
jectory a(i) will be sufficiently close to a eq that Eq. (14. II) can be linearized in the difference 
Aa(t) = a(t) - a eq , to read 

^=Hla>AaW, (4.3) 
with the elements H*ij[a eq ] of the C x C relaxation matrix H*[a eq ] being H*ij[a eq ] = 
(^i[ a ]/ 5o i) a =a^- 

Of course, if we now assume that the fluxes 1Z [a(t)] are indeed given by Eq. (14.21) . then 
we immediately see that the matrix 7i* [a eq ] is given by 

W [a eq ] = -£* [a eq ] oS[a eq ] . (4.4) 
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where the C x C matrix £ [a] is defined as 

and with £ [a eq ] being this matrix evaluated at the equilibrium state a eq . Thus, in the 
neighborhood of the equilibrium state linear irreversible thermodynamics provides a very 
specific structure of the linearized relaxation equation in Eq. (16.51) . namely, 

dAa(t) 



dt 



-r[a>5[a>Aa(t). (4.6) 



The regime in which the relaxation of a(t) is described by this linear relaxation equation 
is referred to as the linear regime around a eq . This is the regime actually discussed by 
Onsager in his fundamental work j^, 4], in which he also discusses the celebrated reciprocity 
relations involving the matrix C* [a eq ] of kinetic coefficients. 



B. Onsager-Machlup theory of thermal fluctuations 

Onsager's connection between the theory of irreversible processes and the theory of spon- 
taneous fluctuations was made by postulating that "the decay of a system from a given 
non-equilibrium state produced by a spontaneous fluctuation Sa(t) = a(£) — a 69 obeys, on 
the average, the (empirical) law for the decay from the same state back to equilibrium, 
when it has been produced by a constraint which is then suddenly removed" This pos- 
tulate, referred to as the Onsager's regression hypothesis, along with the assumption that 
the statistical properties of the fluctuations can be modeled as an Ornstein-Uhlenbeck pro- 
cess, led Onsager, together with his Ph.D. student S. Machlup, to develop their theory of 



time-dependent fluctuations 



l2|. 



Such theory essentially consists of postulating that the dynamics of 5a(t) is described by 
the stochastic version of Eq. (14.61) . namely, 

dSa(t) 



dt 



-L*[a eq ] o £ [a eq ] o 5a(t) + f(t), (4.7) 



where the stochastic vector f(t) is assumed to be a stationary, Gaussian, and purely random 
(or "5-correlated" ) noise, uncorrelated with the initial value Sa(0) of the fluctuations and 
with zero mean and time-correlation function given by < f(t)f T (t') >= j8(t — t'). These 
assumptions, however, are equivalent to saying that 5a(t) is an Ornstein-Uhlenbeck process, 
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as indicated in subsection IIII Al (see Eq. (I3.4p ). Hence, automatically all the mathematical 
properties summarized in section HTT1 apply for the process 5a(t). These include the time 
evolution of the conditional mean value a°(t) and the covariance cr(t), governed by the 
deterministic equations (13. 7h and (13.91) . whose solutions in Eqs. ( \3.6h and (13.101) interpolate 
between the given initial state (a , a ) and the final stationary state (a ss ,a ss ) (which in 
this case corresponds to the equilibrium state (a e<? , <J eq )). These properties, as well as the 
time-dependent correlation function C(t+r, t) in Eq. (13.191) . are thus written in terms of the 
relaxation matrix 7i*[a eq ] or, taking into account Eq. (16. 6p . in terms of the kinetic matrix 
£*[a ei? ] and the thermodynamic matrix £ [a eq \. 

Summarized in the manner we have done in this and in the previous subsection, Onsager's 
theory exhibits a given logical argumentation, in which the foundation of the theory seems 
to be provided by a the phenomenological laws of linear irreversible thermodynamics. In this 
line of reasoning the fluctuations and their description enter as a secondary topic, whose 
description needs just a few auxiliary concepts, including that of an Ornstein-Uhlenbeck 
process. For our purpose of generalizing Onsager's theory, however, an alternative argumen- 
tation is not only conceptually more economical but also easier to extend or generalize. 

This alternative argumentation derives from only two fundamental postulates. The first 
is the assumption that the thermal fluctuations Sa(t) = a(t) — a eq around the equilibrium 
state constitute an Ornstein-Uhlenbeck stochastic process. The second postulate is that 
the first two moments, a ss and a ss , of the stationary distribution Wi[a] of this Ornstein- 



Uhlenbeck process, coincide with 
probability distribution function [62] 



;he first two moments a eq and <j eq of the exact equilibrium 



W eq [a] = eWJ-sia-j;/^ (4. 8 ) 

in which S = 5 [a] is the entropy of the closed system. From these two postulates all the 
results of Onsager's theory follow. 



To see this let us first recall 
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641 ] that from the exact Boltzmann-Planck equilibrium 
distribution in Eq. (14.81) one can derive the value of all its moments, given the fundamental 
thermodynamic relation S = S[a]. In particular one can determine a eq and a eq . Thus, since 
a eq corresponds to a maximum of the entropy, it is determined by the well known and widely 
used first equilibrium condition, 

F[a eq ] = F eq , (4.9) 
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whereas a eq is determined by the less well known but equally important second equilibrium 
condition, 

a eq o£\tfi} = I, (4.10) 

where / is the CxC identity matrix and with the matrix £[a] defined in Eq. (16. 2p . 

The exact condition for the covariance in Eq. (I4.10p can now be used in the mathematical 
result for the structure of the relaxation matrix 7i in Eq. (I3.14p . which then becomes 
identical to the expression for 7i* in Eq. (16.61) . In this manner we identify the matrix L of 
Eq. (13.141) with the matrix £*[a e9 ] of Onsager's kinetic coefficients, and conclude that in the 
present argumentation, the structure of the linear relaxation equations in Eq. (14 .6p can be 
viewed as a consequence of the two fundamental postulates above, and not as a consequence 
of the phenomenological laws of linear irreversible thermodynamics. Instead, we can actually 
invert the conventional argument, and claim that the universal phenomenological validity of 
linear irreversible thermodynamics provides the empirical validation of the suitability of the 
Ornstein-Uhlenbeck model to describe the thermal fluctuations near the equilibrium state. 
An important challenge is now to explore possible manners to generalize Onsager's approach 
to regimes not included in its original formulation. Our proposal to do that is explained in 
the following section. 

V. NON-EQUILIBRIUM EXTENSION OF ONSAGER'S THEORY 

In contrast with the discussion of the previous section, in which we had in mind the sim- 
pler scenario illustrated in Fig. 1(a), let us now consider the schematic scenario illustrated 
by Fig. 1(b), in which we have multiple stationary states corresponding to multiple local 
entropy maxima, and hence, to multiple meta-stable thermodynamic equilibrium states. We 
first consider the application of Onsager's theory to permanently meta-stable equilibrium 
states, then to a ("quasi-static") sequence of meta-stable states in each of which the system 
has sufficient time to equilibrate, and finally to a sequence of meta-stable states in which 
the system does not have sufficient time to equilibrate because the corresponding relax- 
ation times are longer than the time spent by the system in the basin of attraction of the 
instantaneously meta-stable states in the sequence. 
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A. Static entropy landscape and permanent meta-stable states 



Let us assume that the entropy landscape remains static and that, if the system was 
prepared in an initial state a(0) = a with <r(0) = 0, and with a contained in the basin 
of attraction T(a ss ) of the stationary state a ss , the mean trajectory a°(t) will be confined 
to that basin and will relax eventually to the state a ss , in which it will remain indefinitely. 
We then postulate that under these conditions the laws that govern the time evolution of 
the system in this relaxation process are totally indistinguishable from those that govern 
the relaxation towards the most stable equilibrium state a e<? when the system is prepared 
in the basin of attraction T(a e9 ) of this most stable state. Since these laws are contained 
in Onsager's theory, then our postulate is equivalent to the assumption that all what we 
said in the previous section in the context of the state a e<? extends over to all the other 
(meta-stable) equilibrium states & ss without change, except for the replacement of the label 
"eg" by the label "ss". 

For example, the conditional probability density P(a° | a, t) that the system is in the 
macroscopic state a at time t, given that it was prepared in the initial state a e T(a ss ) at 
time t = is given by the Gaussian distribution 

P(a° | a,t) = [(27r) c deta(t)] e -« a - 50 W) to ^ 1 W°( a - 50 ^/ 2 , (5.1) 
where the mean value a°(t) satisfies Eq. (14.61) . which we re- write as 

=f± = -C*[a ss ] o £ [a ss ] o [a°(t) - a ss ], (5.2) 
whose solution is given by 

a°(t) = a ss + e -<c>"K[a"]i o [a - a ss ], (5.3) 

with £ [& ss ] defined as 



v dctj J a=a3S k B \ daidaj , a=aSS 
Similarly, the covariance a(t) satisfies Eq. (13.91) which, using Eqs. (13.111) and (I6.6p . can 
be written as 
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= -C*[a ss ] o S [a ss ] o [a(t) - a ss ] - k B [a{t) - cx ss ] o £ [a ss ] o £*t[ a ss ]. (5.5) 
The corresponding solution, according to Eq. A3. lOj) . is 

a{t) = a ss + e -^[a"]o£[a"]t [(J o _ a s S] Q e (-c^)o S [^))h ^ (5 6) 

which interpolates the covariance a(t) between its initial value cr(0) = a and its stationary 
value a ss given, according to Eq. (I4.10p . by the result 

a ss = (£[a ss ]) _1 . (5.7) 

Finally, according to Eqs. (13.181) and (13.191) . with the matrix 7i* given by Eq. (16.61) . the 
relaxation equation of the two-time correlation matrix C(t + t, t) =< Sa.(t + T)8aJ(t) > is 

dC{t + T ' t} = -£*[a ss ] o S [a M ] o C(t + r, t), (5.8) 
ar 

with initial condition C(t + 0, t) = cr(t), whose solution is 

C(t + r,t) = e - c "^ o£ ^ T o a(t) (for r > 0). (5.9) 

In this manner, if we knew the full thermodynamic landscape described by the fundamen- 
tal thermodynamic relation S = 5 [a], as well as the state dependence of the matrix £*[a] 
of Onsager kinetic coefficients, the results summarized in this subsection would constitute 
the full and exact solution of the Onsager- Machlup model of the irreversible relaxation of a 
closed system towards any stationary state a ss , regardless if this is the most stable or any 
of the meta-stable equilibrium states of the isolated system, but provided that the system 
was prepared in an initial state a that lies inside the basin of attraction of this particular 
stationary state, i.e., that a e T(a ss ). 



B. Quasi-static relaxation to the most stable equilibrium state a eq 

In the previous discussion we assumed that the entropy landscape in Fig. 1(b) remains 
constant in time so that if the system was prepared in the basin of attraction of a given 
stationary state, it will be trapped in this basin, eventually relaxing to the corresponding 
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stationary state a ss , to remain indefinitely in that (stable or meta-stable) state. In real- 
ity, it is impossible that such situation can be sustained indefinitely, since there are many 
possible mechanisms that allow the trajectory a°(t) to explore states outside the basin of 
attraction in which the system was initially prepared. If such excursion occurs to the basin 
of attraction of a stationary state with a still larger local entropy maximum, the second 
law of thermodynamics dictates that the system will proceed to the stationary state of the 
new basin of attraction. If this process occurs repeatedly, the system will eventually reach 
the most stable equilibrium state a eg , thus restoring to some extent the irreversible process 
described in the simpler conditions illustrated in Fig. 1(a), but with a probably much slower 
dynamics because the transitions from basin to basin may involve activation barriers that 
must await for adequate spontaneous fluctuations. 

The irreversible process thus envisioned may in some sense be similar to the concept of 
quasi-static process in classical thermodynamics. Thus, consider a system prepared in an 
initial equilibrium state a/ = a(0) that evolves towards a final equilibrium state a^ through 
a sequence of intermediate equilibrium states a a = a°(t a ) at time t a , with a = 1,2,..., 
aided by the sequential change of the constraints imposed on the system, from those that 
equilibrate the initial state a/, to those that equilibrate the state &i = a°(ti), to those that 
equilibrate the state a2 = a°(t 2 ),..., etc. The time evolution of the covariance o~(t) along this 
process is then described by the interpolating expression in Eq. ( I5.6P which, for the time t 
in the interval t a < t < t a+ i, reads 

a(t) = E-\ x + e C-^ + i<*-+0(*-*-) o [a(t a ) - ° e (-^ + i<^ +1 ) f (*-*<*) 

(t a <t< t a+1 ) (5.10) 

where £ a+1 = £ [a°(t a+1 )] and C* a+1 = C* [a°(t Q+1 )]. 

The main idealized property of a quasi-static process is, of course, that in going from any 
intermediate equilibrium state a Q to the next equilibrium state a Q+ i the intervals | t a+ \ — t a | 
are sufficiently large to allow the system to equilibrate to the new state. This requires these 
intervals to be much larger than any of the relaxation times of the system, i.e., than any of 
the inverse eigenvalues of the relaxation matrix {C* a+l o £ a+ i). The immediate consequence 
of this assumption is that at the end t a+ i of that interval, the covariance matrix a(t) has 
attained its equilibrium value o~(t a+1 ) = £~ +x - Of course, if the system relaxes infinitely fast, 
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we can take these intervals to be arbitrarily small, and then the covariance matrix attains 
its "local equilibrium" value, 

a(t) = {S[sP(t)]y 1 (5.11) 

at all times. A quasi-static process is thus characterized by the assumption that this relation, 
referred to as the local equilibrium approximation, is instantly valid at all times. Thus, for 
example, the two-time correlation matrix C(t+T, t) would be given within this approximation 
and according to Eq. (15. 9p . by 

C(t + r,t) = e -£W)K[s°(t)]T o (£ [a°(i)]) _1 (for r > 0). (5.12) 

In summary, a quasi-static process is characterized by an instantaneously equilibrated 
probability distribution function, i.e., by Eq. (15. ip with a(t) given by its local equilibrium 
value above, so that 

P(a° | a,t) = [(27r) c /det£ [a°(t)]] e H(a-* <'>H*- tf (*»1/ 2 (513) 

Unfortunately, real systems relax through real irreversible processes that have little to do 
with this idealized concept. Instead, the relaxation times assumed negligible in this limit are 
surely finite and, under many circumstances, considerably long. Furthermore, for experimen- 
tal or other reasons, one may want to consider intervals | t a+ i —t a | of much shorter duration 
than the system's relaxation times. Hence, we have in general that a(t) ^ (S [a(£)]) _1 , and 
an important fundamental task is to describe and quantify the deviations from the idealized 
limiting behavior involved in a quasi-static process. The simplest approach to achieve this 
is also based in Eq. (15. 10j) above, as explained below. 

C. Out-of-equilibrium relaxation 

Let us now go back to our closed system, which was prepared in an initial state a that 
is not the global equilibrium state a 69 , but which relaxes irreversibly towards this state of 
maximum global entropy. The statistical description of the macroscopic state of this system 
is provided, as indicated before, by the conditional probability density P(a° | a, t) that 
the system is in the macroscopic state a at time t, given that it was in the initial state 
a at time t — 0. Assuming that the thermal fluctuations <5a(t) = a(i) — a°(t) can be 
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modeled locally in time as an Ornstein-Uhlenbeck process, the Gaussian property of this 
stochastic process approximates P(a° | a, t) by the Gaussian distribution P(a° | a, t) = 
[(2n) c det a(t)]e~ [(a ~ a0( * ))toCT ~ 1( ' )o(a ~ 5 ° ( * )1//2 . Thus, the time evolution of the macroscopic state 
is determined, within this simplified model, by the time evolution of the mean value a°(t) 
and of the covariance <r(t). 

The mean value a°(t) is supposed to solve an equation of the general form in Eq. (14.11) . 
which can also be written as a°(t + At) = a°(t) + 7£[a(t)]Ai + O(At) 2 . Thus, we may 
represent the solution of this equation by a discrete sequence a a = a°(t a ) of the mean value 
at the times t (= 0),ii, ...,t a , ...,too(= oo), generated by the recurrence relation 

a Q+ i = a a + 7Z[a a ]At a , (5.14) 

with the time intervals At a = t a+ i — t a chosen short enough for this linear approximation 
to be valid. 

At this point, we assume that the states & Q of this sequence correspond to instantaneous 
local entropy maxima, and hence, that within the interval (t Q+ i —t a ) the system remains in 
the basin of attraction of the stationary state a a+1 . We also assume that we can approximate 
the corresponding fluctuations 5a(t) = a(t)— a a+ i by an Ornstein-Uhlenbeck process, so that 
in this time interval the time-evolution of the covariance matrix a(t) is described precisely by 
Eq. (|5.10p . This is the expression employed in the previous section to describe quasi-static 
processes. The main difference is that now we are interested in the opposite limit, in which 
the time intervals At a = t a+ i — t a are very short. Thus, we only need to consider Eq. (I5.10p 
in its linear approximation in (t — t a ). The resulting expression leads, with t = t a+ i, to the 
following recurrence relation for the sequence a a = cr(t a ) 



£*a+l ° £a+l O tT a + <T a O S a+1 O C*J +1 (At a ) + £* a+1 + £* T + 1 ( 5 - 15 ) 



where, as before, £ a = £ [a°(t a )] and C* a = C* \aP(t a )\. This, however, is nothing but the 
discrete version of the following differential equation for a(t) 



^ = -£* [a°(t)] o £ [s(t)] o a(t) - a(t) o £ [a°(t)] o C^{t)} + (£*[a°(t)] + C^[a°(t)]) . 

(5.16) 



25 



In contrast with Eq. (15.51) . which has the simple analytic solution of Eq. (15.61) . no analytic 
solution exists for this equation, and hence, the numerical calculation of a(t) may be based 
on the recursion relation in Eq. (I5.15p . The difference between the resulting time dependent 
covariance a{t) and its local equilibrium value {£ [a°(t)]) _1 measures the departure of the 
actual irreversible process from the idealized quasi-static process discussed in the previous 
subsection. 



D. Locally stationary two-time correlation function 

Another important measurable property is the two-time correlation matrix C(t',t) = 



5&(t')5a}(t), where now <5a(t) = a(t) — a°(t) refers to the fluctuations around the time- 
evolving mean value a(t). Thus defined, the fluctuations 5a(t) do not in general constitute 
a stationary stochastic process. However, just like we did in the derivation of the recurrence 
relation in Eq. (I5.15p . here we also assume the local stationarity approximation, i.e., we 
assume that in the interval t a < t < t a+ \, the mean value sf(t a ) and the covariance a(t a ) 
can be considered approximately constant, a°(t) « a°(t a ) and o~(t) ~ cr(t a ), so that the 
thermal fluctuations 5&(t a + r) = a(t a + r) — a°(t a ) around the momentarily stationary 
value 3°(t a ) can be described approximately as an Ornstein-Uhlenbeck process. 

Clearly, the physical notion behind this assumption is that both, a°(t a ) and a(t a ), are 
macroscopic variables that relax to their equilibrium values within rather slow macroscopic 
relaxation times described by the time coordinate t a , whereas the thermal fluctuations 
5a(t a + r) reflect much more local and faster microscopic events, whose correlations de- 
cay within microscopic correlation times described by the time r. These faster events would 
be averaged out when observing only a°(t a ) and a(t a ), and will only "renormalize" the slow 
decay of these properties. Their direct measurement, however, is possible through the mea- 



surement of the r-dependence of the correlation function C(t a + r, t a ) = <5a(t a + r)<5at(t Q ). 

From Eq. (13.171) we may write the stochastic equation for the locally stationary fluctua- 
tions 5&(t a + r) as 

d6a<K \ a + T) = H{t a ) o 5a(t a + r) + f{t a + r) (5.17) 
dr 

with the relaxation matrix 7i{t a ) given, according to the general stationary condition in Eq. 
(EH, by 

H(t a ) = -£(t a )oo--\t a ) (5.18) 
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Multiplying Eq. (15.171) on the right by the transpose of 5a(t a + 0), and averaging over the 
realizations of f(t a + r) and over the "initial" values 5a(t a + 0), we derive the following 



relaxation equation for C{t a + r, t a ) = 5a(t a + r)5at(t c 



dC(t(X ~\~ 7"j /-Q; 



i = -£(* a ) o a- 1 ^) o C(t Q + r, y (5.19) 



whose solution is 

C{t a + r, * Q ) = e - £ (*«)° <T " 1 ^) T o a(t a ). (5.20) 

This expression for C(t a + r,t a ) should be compared with the result in Eq. (15.91) . which 
was derived from Eq. (13.191) with the relaxation matrix 7i given, however, by Eq. (16.61) . 
whose validity is restricted to the linear regime of the stationary state a(t a ). Notice also 
that in the local equilibrium approximation, a{t) = (S [a(£)]) _1 , the result in Eq. (15. 20j) 
coincides with the quasi-static expression in Eq. (15.121) . 



E. Generalized Ornstein-Uhlenbeck processes 

As discussed above, the r-dependence of the thermal fluctuations 5a(t a + r) describes 
faster microscopic events, whose correlations decay within microscopic correlation times. 
When viewed with this temporal resolution, however, it is mandatory to revise the assump- 
tion that the random term f(t a + r) in Eq. (15 . 1 71) can be approximated by a 5-correlated 
noise, an assumption that allowed the dynamics of the thermal fluctuations to be described 
by an Ornstein-Uhlenbeck process. Relaxing this fundamental assumption requires the defi- 
nition of a more general mathematical model of a stationary stochastic process which is not 
required to be necessarily Markovian. 

Such mathematical model, which we will refer to as a generalized Ornstein-Uhlenbeck 



process, was discussed in Ref. [53|] and allows us, in the present context, to describe 5a(t a + 
r) in terms of the most general linear stochastic differential equation with additive noise, 
namely, 

d5a(t a + t) 
dr 

In this model, the stochastic vector f(t a + r) is only assumed to be necessarily stationary 
but not necessarily Gaussian nor white; it is still assumed uncorrelated with the initial value 
5a(t a + 0) and to have zero mean, f(t a + r) = 0. 



f dr'n(t a ;r -r')o6a(t a + r') + f(t a + T). (5.21) 
Jo 
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According to the theorem of stationarity 53j, |54|, the stationarity condition is mathemat 



ically equivalent to the condition that Eq. (I5.2ip conforms to a very strict and rigid format, 
namely, 

d5a{t« + t) = _ w ( ta ) [ .( ta )]-i 0( j a ( ta+r )_y dT'L(t a , r-r')o[(7(t a )]- 1 o5a(t a +r / )+f(ta+r) 

(5.22) 

where u(t a ) is an antisymmetric matrix, u(t a ) = —u)'(t a ), the memory function L(t a ;r) 
satisfies the following fluctuation-dissipation relation 

L{t a - r) = L\t a - -t) =< f{t a + r)f \t a + 0) >; (5.23) 



the matrix o~(t a ) = 8a.(t a + 0)5al(t a + 0) is the covariance of the probability distribution 
function of the initial value 5a.(t a + 0) of the fluctuations. In other words, the non-Markovian 
relaxation matrix TC(t a ; r) of Eq. (15.211) is not arbitrary, but must have the rigid format 
leading to Eq. (I5.22p . namely, 

H(t a ; t) = - [25(r)uj(t a ) + L(t a , r)} o [a(t a )}- 1 . (5.24) 

The generalized Ornstein-Uhlenbeck model defined by Eq. (15.221) may be used to describe 
the fluctuations around any stationary state, including meta-stable and absolutely stable 
thermodynamic equilibrium states. Under stationary conditions the label t a is, of course, 



unnecessary and Eq. (15.221) may be recognized in a statistical mechanical context 51[ as the 
generalized Langevin equation (GLE). In that context, however, the term GLE is associated 
with the stochastic equation formally derived from a N-particle microscopic (Newtonian 
or Brownian) dynamic description by means of projection operator techniques to describe 



the time-dependent thermal fluctuations of systems in thermodynamic equilibrium [51]. 
Indeed, such an equation conforms exactly to the format described by Eq. (I5.22p . It is 
important to insist, however, that this format has a purely mathematical origin, imposed 
solely by the stationarity condition, and is certainly NOT a consequence of the formal 
possibility of deriving it from an underlying microscopic level of description. In fact, it is 
this mathematical structure of the GLE, and the "selection rules" imposed by the symmetry 
properties of the matrices u ( t n ) and L(t a ;r) (along with other selection rules imposed by 



other possible symmetries 53(), what allows a fruitful use of the rigid format of this equation 



to describe complex dynamic phenomena in a rather simple manner, with virtually complete 
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independence of the detailed N-particle microscopic dynamics underlying the time-evolution 
of the fluctuations Sa(t). 

Multiplying Eq. (15.221) on the right by 5&^(t a ), and taking the corresponding average, we 



have that the two-time correlation function C(t a + r, t a ) = 5a(t a + r)5at(t a ), that we shall 
now denote as Ct a ( T ), satisfies the following relaxation equation 

dC ta (r) 



dr 



-u(t a ) o [a^)]- 1 o C ta (r) - [ T dr'L(t a , r - r') o [a^)]" 1 o C ta (r') (5.25) 

Jo 



whose solution may be written, in terms of the Laplace transforms C ta (z) and L(t a ;z) of 
C ta (r) and L(t a ; r), as 

C ta (z) = {zl+ [u(t a ) + L(t a , z)\ o [ait^y 1 o a(t a ). (5.26) 

Let us notice that the results in this subsection reduce to the results of the previous 
subsection in the Markovian limit, defined by the condition L(t a ; z) w L(t a ; z — 0) = L°(t a ) 
or, equivalently, by the condition L(t a ; r) w 26(r)L°(t a ), with L°(t a ) = J °° dtL(t a ; r). Thus, 
it is not difficult to see that Eqs. (I5.22p . (I5.25p . and (I5.26P above become, respectively, Eqs. ( 
I5.17p . (I5.19p . and (15.201) of the previous section, and that Eq. (I5.24p for TC(t a ; r) corresponds, 
after integrating over r, to Eq. (15.181) for 7i{t a ) upon the identification of the matrix C(t a ) 
with 

C(t a )=u(t a ) + L°(t a ). (5.27) 

Finally, let us postulate a "correspondence principle" that guarantees that in the vicinity 
of the stable thermodynamic equilibrium state a 69 , the present nonlinear theory reduces to 
Onsager's original linear theory. For this, we assume that the matrices w(t a ) and L(t a ; r) 
depend on the relaxation time t a only through a(t a ), so that uj(t Q ) = u[a(t a )} and L(t a ; r) = 
L[t; a(toi)]. Then, Eq. ( 15.271) may be rewritten as 

POO 

C(t a ) = uj[a(t a )] + / dTL[r;a(t a )}. (5.28) 
Jo 

We then postulate that C(t a ) must coincide, when the system has fully relaxed to the equi- 
librium state a 69 , with the phenomenological matrix £*[a e? ] of Onsager's kinetic coefficients 
involved in the linear laws of Eq. (14. 6p . i.e., that 

POO 

£*l& eq ] =u[a eq ]+ / drL[r;si eq }. (5.29) 
Jo 
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VI. GENERAL DESCRIPTION OF OUT-OF-EQUILIBRIUM RELAXATION 



In this section we recapitulate the discussion of the previous section as a final pro- 
posal for the general canonical description of the nonlinear non-equilibrium relaxation of 
macroscopic systems. This scheme is summarized by the time-evolution equations for the 
mean value a(t), for the covariance cr(t), and for the time- dependent correlation function 



C t {r) = 5a{t + T)Sa^(t). 

Let us first mention that an essential piece of information that must be provided exter- 
nally to this canonical theory is the fundamental thermodynamic relation S = S[a\. From 
this relation one is supposed to determine the thermodynamic equations of state, i.e., the 
functional dependence on a of the intensive variables 

Fi[a] = k B \dS[ a ]/d ai ), (6.1) 

as well as the thermodynamic matrix £ [a], defined as 

\ dcij ) k B i^daida^) ' ^ 

A stable thermodynamic equilibrium state is determined by the equilibrium condition for 
& eq in Eq. (14.91) . namely, 

F[a e "] = F eq , (6.3) 

and by the condition in Eq. (14.101) that determines the equilibrium value a eq of the covari- 
ance, namely, 

a eq o£[ a eq ] = I, (6.4) 

In the neighborhood of a thermodynamic equilibrium state, the generally nonlinear re- 
laxation equation for the conditional mean value a°(t) can be linearized in the difference 
Aa(t) = a°(t) - a eg , to read 

^=Hla>AaW, (6.5) 
with the matrix 7i* [& eq ] given by 

H*[a eq ] = -C*[a eq ]o£[ a eq ], (6.6) 

with £*[a eg ] being the matrix of Onsager's kinetic coefficients. 
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Outside this linear regime, however, the mean value a(t) is governed by the phenomenolog- 
ical transport equation for the corresponding macroscopic state variables, written in general, 
according to Eq. (14. ip . as 

d ^ = nw)\- (6.7) 

The constitutive relation, which determines the generally nonlinear dependence of the vector 
1Z [a(t)] on the state vector a(t), is in general unknown and is regarded as another external in- 
put of the present theory. In the vicinity of a thermodynamic equilibrium state a eq , however, 
this non-linear dependence must reduce to the linear laws of irreversible thermodynamics. 

Assuming that the matrices C*[a] and £ [a] are defined at any accessible state a, the 
time-evolution of the covariance matrix a(t) is governed by Eq. (15.161) . i.e., 

^ = -£*[a°(f)] o £ [a(f)] o a(t) - a{t) o £ [a°(t)] o £*[a°(*)] + (£*[a°(t)] + £*[a°(*)]) , 

(6.8) 

where a°(i) is the solution of the nonlinear equation (16.71) . This equation for a(t) has no 
explicit solution, and hence, must be solved simultaneously with Eq. (16 .7p for a(i) using, 
for example, the recursion relations in Eqs. (15.141) and (15.151) . 

The dynamics of the locally stationary fluctuations Sa.(t + r) = a(i + r) — a°(t) around 
the conditional mean value a°(i) are described by Eq. ( 15.22ft . which reads 



d6a ^ + T) = -wja^oa^ttlo^ + r)- f dr'L[r- r'; a(t)] o^ 1 ^) o5a(t + r') +f(t + r), 
c* 7 " Jo 

(6.9) 

with < f(t + r) >= and < f(t + r)^ (t + r') >= L[ t — t'\ a(t)]. From this equation one 
derives the time-evolution equation of the non- stationary time correlation function Ct(r) =< 
5a(t + r)Sa^(t) >, which reads 

OC t (r) 



Or 



-w[a(t)] oa" 1 ^) oC t (r) - / rfr'L[r - r'; a(t)] o ^(t) o C t (r'), (6.10) 

Jo 



whose initial condition is C t {j = 0) = cr(t). This equation describes the decay of the 
correlation function C t {r) with the "microscopic" correlation time r, after the system has 
evolved during a "macroscopic" relaxation time t from an initial state described by a = 
a(t = 0) and a = a(t = 0), to the "current" state described by a(t) and <r(t). 
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Finally, we postulate a connection between the phenomenological matrix £*[a ei3 ] of On- 
sager's kinetic coefficients involved in the linear laws of Eqs. (16.51) and (16. 6p . and the matri- 
ces a; [a] and L[t; a] that describe, according to Eq. (I6.9p . the dynamics of the fluctuations 
5a(t + r). This connection, established in Eq. (I5.29P above, reads in general 

POO 

C*[a}=cu[a] + drL[r;a], (6.11) 
Jo 

and constitutes a correspondence principle which ensures that the theory just summarized 
contains the conventional Onsager's theory as a particular case in the vicinity of a thermo- 
dynamic equilibrium state. 

In this manner, for given initial conditions a(t = 0) = a , a(t = 0) = cr°, and Ct{r = 
0) = <r(t), Eqs. (IP) . (RTSj) . and flBTTUD . together with the relationship in Eq. (I6TTD . 
would constitute a closed system of equations if two fundamental pieces of information were 
available. The first is the fundamental thermodynamic relation S = S[a], from which the 
state-dependence of F[a] = (dS[a]/da) and £ [a] = — (<9F[a]/<9a)/&B could be derived. The 
second refers to the conservative and dissipative kinetic matrices, o;[a(i)] and L[r;a(t)], 
entering in Eq. (16.101) . These two fundamental pieces of information must be provided 
externally to the general format above, and in many cases, their investigation constitutes 
a relevant problem by itself. However, for a given specific physical context, the format 
just described may guide us in the construction of the specific models and approximations 
that best suit the description of a particular relaxation phenomenon. For example, the 
application of the generalized Onsager's theory just summarized, to the specific context of 



52] . In order to illustrate 



the dynamics of colloidal dispersions, was spelled out in Ref. 
the concrete use of the general theory in a concrete physical condition, in what follows we 
present a brief review of some aspects of such application. 



VII. APPLICATION TO COLLOID DYNAMICS 



In this section we apply the general concepts above, to the specific problem of the dif- 
fusive relaxation of the local concentration of particles in a colloidal dispersion without 
hydrodynamic interactions. Thus, let us consider a dispersion of N colloidal particles of 
mass M in a volume V which, in the absence of external fields, has a uniform bulk number 
concentration ris = N/V. In the presence of a conservative static external field that exerts 
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a force F ex '*(r) = — V\l/(r) on one particle located at position r, the mean local concen- 
tration profile of colloidal particles, n(r,t), will evolve in time from some initial condition 
n(r,t = 0) = n° (r) towards its stable thermodynamic equilibrium value n eg (r). The initial 
profile n° (r) is, of course, arbitrary, whereas the final equilibrium profile n eq (r) is univocally 
dictated by the external and internal constraints on the system according to the second law 
of thermodynamics. In practice, the external constraints are represented by the potential 
^(r) of the external forces, while the internal constraints originate in the intermolecular 
interactions, represented by a pair potential u(r, r'). 

To simplify the correspondence with the general theory of the previous section, let us 
imagine that we divide the volume V occupied by the colloidal dispersion in C cells of equal 
volume AV^ fixed in space. We then describe the macroscopic state a(t) of this system in 
terms of the variables eij(t) = n^t) = Ni(t)/AV, where Ni(t) is the number of colloidal 
particles in cell % at time t. We shall employ the results of the previous section, however, in 
the understanding that the continuum limit, C — > oo and AV — > 0, has been taken. In this 
limit, the label i (= 1, 2, C) of the component a, is changed to the label r e V denoting 
the spatial position of the center of the cell, and the component a,i(t) becomes the function 
n(r,t), which is the local concentration profile of colloidal particles at time t. 

Let us first discuss the application of the general equilibrium conditions in Eqs. (I6.3P and 
( 16. 4p that determines, in the equilibrium state, the mean value n eq (r) and the covariance ma- 
trix a eq (r, r'), and then we identify the kinetic information by means of a phenomenological 
derivation of the generalized diffusion equation. 



A. Fundamental thermodynamic relation and equilibrium conditions 

The fundamental thermodynamic relation of the present system expresses the functional 
dependence of the entropy S on the local concentration profile n(r), a dependence repre- 
sented by S = S[n]. The intensive variable conjugate of n(r) is — /3{i[r;n], i.e., it is the 
negative of the local electrochemical potential /x[r; n(t)} at position r in units of the thermal 
energy k^T = f3~ x . It is an ordi nary function of r and a functional of the concentration 



profile n(r), written in general as 65] 
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Pufrn] = pfi m [r-n]+^(r) (7.1) 
= Pn\P) + lnn(r) - c[r; n] + 0* (r). 

The first two terms of this definition of // m [r;n], (/3/i*(/3) + lnn(r)), are the ideal gas con- 
tribution to the chemical potential, whereas the term — c[r; n] contains the deviations from 
ideal behavior due to interparticle interactions. 

The equilibrium condition in Eq. ( 16.31) . which determines n eq (r), may then be written as 

P/i[r;n eq ] = /3fi*(P) + lnn eq {r) - c[r; n eq ] + /3tf(r) = (7.2) 

where the constant fi R is the uniform value of the electrochemical potential. This would be 
a closed equation for n eq (r) if we knew the functional dependence of c[r;n] on n(r). The 
simplest example of a proposed approximate functional dependence is the linear functional 
c[r;n] = — J dr'/3u(r,r')n{r'), referred to as the Debye-Hiickel approximation. 

The thermodynamic matrix £ [a] defined in Eq. f l6.2j) can be written in general, using 
Eq. flU}, as 



£[t, r'; n] 



5(r-r')/n(r) - c (2) [r, r'; n], (7.3) 



8n(r') 

with [r, r'; n] = (Sc[r; n\/8n(v')) being the functional derivative of c[r;n] with respect 
to n(r'), referred to as the direct correlation function. On the other hand, the covariance 



matrix cr(r, r') = 5n(r, 0)5n(r', 0) can be written in terms of the total correlation function 



h^(r ,r') as 



o~(r, r') = n(r)5(r — r') + n(r)n(r')/i < - 2 - ) (r, r'). (7.4) 

The matrices S[r, r'; n] and o~(r, r') are not in general related to each other. When evaluated 
at the equilibrium local concentration profile n eq (r), however, they are related by the second 
equilibrium condition in Eq. (16.41) . which in the present context reads 

J dr'a eq (r, r')£[r', r"; n eq ] = 5{r - r"). (7.5) 

One can immediately see that this equation is equivalent to the well-known Ornstein-Zernike 
equation 
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h(r, r') = c(r, r') + / dV'c(r, r")n eq {v")h{r" , r'). (7.6) 

where c(r, r') and h(r,r') are, respectively, the equilibrium value of c^(r, r') and h^ 2 '{r, r'). 
Sometimes approximate chemical equations of state may be expressed as a "closure" relation 
between these two properties, as in the so-called hyper-netted chain (HNC) approximation, 



which writes 



12| 



c(r, r') = -/3u{r, r') + h{r, r') - ln[l + h(r, r')]. (7.7) 

Within this approximation one would have to solve self-consistently Eqs. (17. 2p . (17.61) . and 
(17. 7j) . together with c[r; n eiJ ] = J dr'c(r, r')n eq (r'), for the properties n e9 (r), c[r;n eg ], c(r, r') 
and /i(r, r'), given the pair potential u(r, r') and the potential \I/(r) of the external field. 

B. Non-equilibrium diffusion in colloidal dispersions 



In Ref. 52] the time-evolution equation for the local concentration profile n(r,t) was 
derived by complementing the continuity equation, 

^ = -V-J(r,t), (7.8) 
with a constitutive relation constructed by defining the particles' velocity field u(r, t) by 

u(r,t)=j(r,t)/n(r,t). (7.9) 
Thus, u(r, t) is the velocity of a particle representative of the set of particles in a volume 



dr centered at position r. One may then essentially follow Einstein's proposal [66( that the 
friction force on this representative particle must be equilibrated, on the average, by the 
osmotic force — V/z m [r;n] and by the external force — V\l/(r), both included in — V//[r;n]. 

The friction force is the sum of the friction due to the supporting solvent, — C°u(r, t), and 
the frictional effects of the interactions with the other particles, — (A£)u(r, t). The latter, 
however, may involve in general spatial and temporal nonlocal effects due to the collective 
character of u(r, t), so that the static equilibrium of these forces must actually be written 
in general as 

C°u(r, t) + J dt'J d 3 r'AC(r -r'-t- t')u{r', f) = -V/i[r; n{t)]. (7.10) 
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whose solution for u(r, t) is given by 

u(r, t) = —D° J dt' j d 3 r'b[r -r';t- t']V(3(i[r'- n{t% (7.11) 

with the spatially and temporally non-local collective mobility kernel b[r — r';t], viewed as 
the (r, r') element of the "matrix" b(t), being the solution of a "matrix" equation which in 
Laplace space reads 

b(z)o [J + ACW/C ] = 1 (7.12) 

where the matrix product " o " means spatial convolution, and "J" is Dirac's delta function 
8(r — r'). In this equation, the matrix b(z) = J °° dte~ zt b{t) is the Laplace transform of 
the matrix b(t), and similarly for A((z). The free diffusion coefficient D° is defined as 
D° = k B T/(°. 

Using Eq. (17.111) in j(r, t) = u(r, t)n(r, t), the continuity equation (17. 8p finally leads us 
to the most general diffusion equation, namely, 

dn(r,t) ''' 
dt 

Let us now discuss the use of this equation to describe the relaxation of the macroscopically 
observed mean value n(r, t) and of the covariance a(r, r'; t), as well as of the time- dependent 
correlation function C(t,t') of the fluctuations 5n(r,t) = n(r,t) — n(r,t). 



D°V ■ n(r, t) J dt' J d 3 r' b[r -r';t- t']V [3 fi[r' ; n(t')} . (7. 13) 



C. Irreversible relaxation of n(r,i), er(r, r';i), and C(t,t') 

Following the general format of the relaxation equations presented in section IVT1 (i.e.. Eqs. 
(16. 7p . (16. 8p . and (16. 9p . and (I6.10p ). we assume that the spatial and temporal arguments of 
the mean local concentration n(r, t) describe spatial and temporal variations of macroscopic 
scale, so that, for example, in a quenching process, the variable t is the ageing or waiting 
time. In contrast, the thermal fluctuations 5n(r, t + r) = n(r,t + r) — n(r,t) vary within 
microscopic times denoted by r which may be much shorter than t. In a similar manner, we 
also assume that the spatial variation of n(r, t), described by the spatial argument r, occur in 
much larger spatial scales than the microscopic spatial variations of the thermal fluctuations 
5n(r + x, t + r) indicated in the neighborhood of r by the spatial coordinate x. Thus, 
our central assumption is that the mean value n(r, t) remains approximately uniform and 
stationary while the fluctuations vary microscopically within the finer space and time scales 

36 



indicated by x and r. We express this assumption by describing the macroscopic relaxation 
of n(r,t) by the temporally and spatially local version of Eq. (17.131) . This corresponds to 
approximating the generalized mobility kernel b[r — r';t — t'] of this equation by 

b[r-r';t-t'} = b*{r, t)5{r - r')25{t - t'), (7.14) 

where 

/poo 
rfx J dr 6[x,r;r,t] (7.15) 

with 6[x, r; r, t] = b[(r + x) — r; (t + r) — t]. In this manner we can write the analog of Eq. 
(16. 7p . i.e., the diffusion equation for the mean value n(r, t), which reads 

We may now linearize this equation around the equilibrium profile n eq (r), to get the 
analog of Eq. (16. 5p and (16. 6p . From the resulting linearized equation we can identify the 
"matrix" C*[r, r'; n(t)] of Onsager kinetic coefficients as 

- C*[r,r';n(t)} = D°V ■ n(r, t) b*(r, t) V5(r - r). (7.17) 
and, from Eq. (16.81) . we can write the relaxation equation for cr(r, r'; t) as 



dt 



da(r,r';t) =£) o v . ^ f) 6 *( rjt ) V J dr 2 £ [r, r 2 ; n{t))a{r 2 , r'; t) 

+ D°V ■ n(r', t) b*{r', t)V J dr 2 £[r', r 2 ; n(t)]a(r 2 , r; t) 



(7.18) 



- 2L>°V • n(r, t) 6*(r, t)V5(r - r'). 

Also according to the generalized Onsager scheme, the dynamics of the fluctuations 
8n(r,t + r) = n(r,t + r) — n(r, t) are now described by a stochastic equation with the 
structure of Eq. (16.91) . In our case, this equation is meant to describe the relaxation of the 
fluctuations 5n(r,t + r) in the temporal scale described by the time r, around the mean 
value n(r, t) of the local concentration at position r and time t. The assumption of local 
stationarity means that in the time-scale of r, n(r, t) is to be treated as a constant. Although 
not explicitly contemplated in the format of Eq. (16.91) . but as already indicated above Eq. 
(I7.14p . here we also add the spatial counterpart of this simplifying assumption. Thus, we 
write the fluctuations as Sn(r + x, t + r) = n(r + x, t + r) — n(r, t), where the argument r of 
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n(r, t) refer to the macroscopic resolution of the measured variations of the local equilibrium 
profile, whereas the position vector x adds the possibility of microscopic resolution in the 
description of the thermal fluctuations. Defining the fluctuations as the deviations of the 
microscopic local concentration profile n(r + x, t + r) from the mean value n(r, t) indicates 
that, within the microscopic spatial variations described by the position vector x, n(r, t) 
must be treated constant. 

With this understanding, we can now proceed to identify the elements of Eq. (16. 9p 
corresponding to our problem. In the present case, the corresponding antisymmetric matrix 
u;[a(£)] vanishes due to time-reversal symmetry arguments [53|. We can then write the 
matrix L[r;a(t)] as the non-markovian and spatially non-local Onsager matrix implied by 
the general diffusion equation in Eq. ( 17. 13ft . which must reflect, in addition, that within the 
temporal and spatial resolution of the variables x and r, the local concentration profile n(r, t) 
remains uniform and stationary. These assumptions can be summarized by the following 
stochastic equation for 8n(r + x, t + r) 



d5n(r + x, t + r) 



Or 



:D°n(r,£)V x • / dr' f dxi6[x - xi, r - r'; r, i]V Xl 

Jo J (7.19) 

/ dx 2 o- _1 (xi, x 2 ; t)6n(r + x 2 , t + r') + f(r + x, t + r), 

where the function (j _1 (x, x'; t) is the inverse of the covariance <r(x, x'; t) in the sense that 

J rfx" ^(x, x"; t)a{x", x'; t) = <J(x - x). (7.20) 

The random term f(r + x, t + r) of eq. (I7.19P is assumed to have zero mean and time 
correlation function given by < f(r + x, t + t)v (r + x ; , t + r') >= L[x — x', r — r'; r, i], with 

L[x - x', r; r, t] = D°n(r, t)V x ■ / rfx^fx - Xi, r; r, t] V Xl <5(xi - x'). (7.21) 

Similarly, the analog of Eq. (I6.10p for the time correlation function C t {r) is the relaxation 
equation for C(x — x', r — r'; r, t) =< 5n{r + x, t + r)5n(r + x', t + r') >, namely, 

ag(x r;r,t) =£>0 _^ f) ^ f _ XljT _ /.^^ 

° T Jo J (7.22) 

/ dx 2 cr _1 (xi, x 2 ; t)C(x 2 , r'; r, t). 



38 



In this manner the generalized theory of non-equilibrium diffusion just presented writes 
the relaxation of the mean value n(r,t), of the covariance a(r, r';i), and of the time- 
correlation function C(x, r; r, t), through Eqs. ( 17.161) . (I7.18p . and (17.221) . in terms of the 
generalized mobility 6[x, r; r, t] or, according to Eq. ( 17.121) . in terms of the temporally 
and spatially nonlocal friction function A£[x, r;r, £]. These equations constitute the gen- 
eral framework in which approximations may be introduced to construct a closed system of 
equations for the properties involved. This was in fact the main aim of Ref . [52| , and hence, 
at this point we refer the reader to that reference for further details. 

VIII. CONCLUDING REMARKS 

In this paper we have proposed a generalization of Onsager's theory of the time-dependent 
thermal fluctuations <5a(£) = a(t) — & eq around the equilibrium state a e<? , to the description 
of the thermal fluctuations <5a(t) = a(t) — a(t) around the time- dependent mean value a(t) 
that relaxes irreversibly towards its most stable equilibrium state a. eq as the solution of a 
(generally nonlinear) relaxation equation. The essential results of this generalized theory 
were summarized in section IVfl and consist of the relaxation equations for the covariance 
matrix a(t) = <5a(i)<5a(t) in Eq. (16.81) and for the two-time correlation function Ct(r) = 
5&(t + r)5a(t) in Eq. ( 16.101) . The time t represents the macroscopic relaxation time that 
describes the time-evolution of both, the mean a(t) and the covariance cr(t), whereas the time 
r in C t {j) represents the ("microscopic") correlation time of the fluctuations, as observed 
at the macroscopic time t after the system was prepared at the given initial conditions 
a(i = 0) = a and a(t = 0) = a . 

Eqs. ( 16. 7p . ( 16. 8ft . and ( 16.101) . together with the relationship in Eq. ( 16. lip , would constitute 
a closed system of equations if two fundamental pieces of information were available. The 
first is the fundamental thermodynamic relation S = 5 [a] and the second refers to the 
conservative and dissipative kinetic matrices, u;[a°(£)] and L[t; a(t), cr(t)], entering in Eq. 
( I6.10p . These two fundamental pieces of information must be provided externally, and must 
reflect the specific context of a particular relaxation phenomenon. In order to illustrate the 
use of this generalized canonical theory by means of a particular application, in the previous 
section we described its application to the specific context of the dynamics of colloidal 
dispersions. There we first discussed the main features of the fundamental thermodynamic 
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relation of these systems and then proposed physical arguments leading to a generalized 
diffusion equation. 

These are the two elements that had to be provided externally to the canonical theory. 
With these two inputs discussed, and following the script of the canonical theory, we could 
write the time evolution equations for the mean local concentration profile n(r, t) and for the 
covariance <r(r, r'; t) of the thermal fluctuations <5n(r, t) = n(r,t) — n(r,t) (Eqs. (17.161) and 
(I7.18p . respectively). The time evolution of these two properties, as they relax irreversibly 
from some arbitrary initial values rf(r) and <T°(r,r') to their equilibrium values n eq (r) and 
a eq (r, r'), constitute the fundamental results for the description of the irreversible relaxation 
of the system. 

We then derived the relaxation equation in Eq. (I7.22p for the time-dependent correlation 



function C(x, r; r, t) = 5n{r + x, t + r)<5n(r, t) in terms of the generalized time-dependent 
friction function A£[x, r; r, t\. This is in fact the memory function of a generalized mobil- 
ity function 6*[x, r; r, t] appearing in Eq. ( 17. 221) and, in its Markovian limit, also in Eqs. 
(IT. 161) and (17. 18ft . Thus, in summary, in the previous section we expressed the main proper- 
ties that describe the dynamics of the non-equilibrium relaxation of a colloidal dispersion, 
namely, n(r, t), cr(r, r';t), and C(x, r;r, £), in terms of the memory function A£[x, r; r, t}. 
The corresponding system of equations constitute the non-equilibrium generalization of the 
fundamental equations upon which non-equilibrium generalization of the SCGLE theory 
of the dynamics of equilibrium dispersions can be constructed. In fact, Eq. (17.221) for 
C(x, r; r, t), together with Eq. (17.121) . is the non-equilibrium analog of the memory equa- 
tion for the intermediate scattering function F(k,t) in terms of the so-called irreducible 
memory function A((k,t), as indicated in detail in Ref. 52j. There a proposal is presented 
of the full non-equilibrium generalization of the SCGLE theory of colloid dynamics and of 
its application to dynamic arrest phenomena. 

In the same reference it is pointed out that the same non-equilibrium scheme summarized 
in the previous section contains as particular cases some results and equations that are 
important in various specific contexts. For example, the diffusion equation for n(r, t) in Eq. 
(I7.16p . in the limit in which we ignore the memory effects contained in A((k,t), happens 
to coincide with the central general equation of the recently-developed dynamic density 



functional theory [67, 



681 ]. which has been applied to a variety of systems, including the 



description of the irreversible sedimentation of real and simulated colloidal suspensions [69 ]. 
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It also coincides in certain circumstances with an equation for the irreversible relaxation of 



n(r,t) derived by Tokuyama 70j, |7l|. On the other hand, the relaxation equation for the 



covariance a(r, r';i) in Eq. (17.181) . also in the limit A((k,t) = 0, may also be shown to 



contain, as a particular case, the fundamental equation emp 
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oyed in the classical description 



of the early stages of spinodal decomposition [72, |73|, 

Thus, the general results of the previous section can be used as the basis for the extension 
of the aforementioned theories to the case in which the memory effects due to the direct 
interactions, contained in A((k,t), cannot be ignored. These effects are responsible, for 
example, for the dynamic slowing down of the system and for its eventual transition to 
dynamically arrested conditions. Thus, it is this more general scheme that is expected to 
generate the most original predictions including, for example, the dependence of the glass 
transition scenario on the cooling rate or the description of the ageing of the static structure 
factor and of the intermediate scattering function after a quenching process. Specific and 



more concrete advances in this direction, however, will be reported separately 76]. 
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